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ABSTRACT 

The candidate SXPhe star KIC 11754974 shows a remarkably high number of com- 
bination frequencies in the Fourier amphtude spectrum: 123 of the 166 frequencies 
in our multi-frequency fit are linear combinations of independent modes. Predictable 
patterns in frequency spacings are seen in the Fourier transform of the light curve. We 
present an analysis of 180 d of short-cadence Kepler photometry and of new spectro- 
scopic data for this evolved, late A- type star. We infer from thell50-d, long-cadence 
light curve, and in two different ways, that our target is the primary of a 343-d, 
non-eclipsing binary system. According to both methods, the mass function is simi- 
lar, f{M) = 0.0207 ± 0.0003 Mq. The observed pulsations are modelled extensively, 
using separate, state-of-the-art, time-dependent convection (TDC) and rotating mod- 
els. The models match the observed temperature and low metallicity, finding a mass 
of 1.50-1.56 Mq. The models suggest the whole star is metal-poor, and that the low 
metallicity is not just a surface abundance peculiarity. This is the best frequency anal- 
ysis of an SX Phe star, and the only Kepler 5 Set star to be modelled with both TDC 
and rotating models. 

Key words: Stars: oscillations - stars: variables: delta Scuti - stars: Population II - 
stars: individual (KIC 11754974) - physical data and processes: asteroseismology. 



1 INTRODUCTION 
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The goal of the Keple r mission is to fin d Earth-like planets 
in the habitable zone (JKoch at al.ll2010l ). The methodology 
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inherent in this task is the simultaneous, space-based ob- 
servation of ~150 000 stars in a 115 deg^ field of view for 
transiting events. 

Kepler observations are made in two cadences: most 
observations take place in long cadence (LC) format, where 
photometric readouts are co-added for 29.4 min. Short ca- 
dence (SC) observations have effective integration times of 
58.8 s. These cadences make the data excellent for asteroseis- 
mic investigations. Kepler data are organised into quarters, 
denoted Qn, and SC data are further subdivided into three 
one-month segments, denoted Qn.m. 

The success of the Kepler mission lies in its unprece- 
dented /imag photometric precision. Even for a star as faint 
as KIC 11754974 (Kp = 12.7 mag), the subject of this work, 
the average noise amplitude in the Fourier transform is 
~5 /iinag at low frequency, decreasing gradually to ~ 2 /xmag 
for frequencies above about 70 d~^. 

Kepler photometry not only offers higher precision, but 
also provides nearly-continuous observations. Progress in the 
analysis of 5 Set stars is no longer hampered by the insur- 
mountable time gaps in the data associated with ground- 
based observations. Most stars now have almost continuous 
observations in LC mode from the commissioning run of QO 
right through to the most recent quarter at the time of writ- 
ing (Q13), spanning over 1150 d; the same is true for some 
targets in SC mode. This long time-span is particularly im- 
portant for resolving peaks that are very closely spaced in 
frequency, and for studies of frequen cy or amplitude modu- 
lation, among other reasons (see, e.g.. lMurphvl2012bl '). Com- 
pletely alleviating Nyquist am biguities with Kepler data can 
now a lso be added to that list (jMurphv. Shibahashi fc Kurta 
I2OI3I ). 

The 5 Set stars are located at the junction of the clas- 
sical Cepheid instability strip and the main-sequence, with 
pulsation periods between 18 min and 8hr. They are driven 
by the opacity mechanism in the Hell partial ionisation 
zone. SX Phe stars are the Population II counterparts of the 
5 Set stars, charac terized by low metallicity and high tan- 
gential velocities (|Nemec fc Mated 1 1990: Bal ona fc Nemed 
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I2OI2 V Within the last decade, observations of SXPhe stars 
have revealed that the y are not necessarily intrinsically high- 
ampl itude pulsators (lOlech et al.l l2005l : iBalona fc Nemed 
[2013), but show amplitude distributions similar to the Pop- 
ulation I 5 Set stars, viz., having peak pulsatio n amplitudes 
of a few mmag (see lUvtterhoeven et al.ll201ll for a statisti- 
cal study) . Traditionally, 5 Set stars with amplitudes above 
0.3 mag were termed HADS (High Amplitude Delta Scuti) 
stars, and were found to pulsate mostly in one or two ra- 
dial modes. Their high amplitudes are attributed to their 
evolved stages (e.g. [Petersen fc Christensen-Dalsgaardlll996l . 
who describe HADS stars as being in the immediate post- 
MS stage). KIC 11754974 has a maximum peak-to-peak light 
variation of only 0.24 mag; base d on tangential ve l ocitie s was 
identified as an SXPhe star bv lBalona fc Nemed (|201^) and 
was the highest-amplitude Kepler SXPhe star they identi- 
fied. In this work, we find KIC 11754974 to have properties 
in common with both HADS and SXPhe stars. 
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Figure 1. Temporal changes in the frequency of the four 
strongest modes excited in KIC 11754974, /i, /2, /a and /4, top 
to bottom. The values derived from LC and SC data are plotted 
with filled and open symbols, respectively. The sinusoid with a 
period of 343 d fitted to the /i data is also shown. For clarity, an 
offset has been added to the /2, /a and /4 data; the sinusoid has 
also been shifted accordingly. The right-hand side ordinate shows 
radial velocities calculated from the Doppler effect. 



2 BINARITY OF KIC 11754974 

2.1 Doppler-shifted frequencies 

Our target falls on Kepler's failed Module 3, meaning it is 
only observed for three-quarters of each year. For our fre- 
quency analysis we utilized the higher sampling of the SC 
data, the longest uninterrupted run of which is Q6 and Q7, 
covering about 180 d. In the course of analysing the com- 
bined Q6-|~Q7 dataset it appeared that in Fourier spectra 
a significant residual signal occurs near the strongest terms 
after their subtraction. This indicated frequency and/or am- 
plitude changes which prompted us to carry out the analysis 
separately for shorter datasets. Both LC and SC data were 
used in this analysis; they were split into 30-50 d subsets. For 
each subset, the frequency set containing the 72 strongest 
terms (20 independent and 52 harmonics and combinations) 
was fitted by means of non- linear least-squares. The result- 
ing frequencies of the strongest modes in different subsets 
were then compared. The result is shown in Fig.[T] where we 
plot the fractional frequency change of the four frequencies 
with the highest amplitudes in a given time subset, relative 
to their amplitudes in the arbitrarily chosen reference sub- 
set, here the SC Q6 data. As can be seen, the frequencies of 
all four of the strongest modes behave in the same manner 
and can be described by a sinusoid with a period of about 
343 d. 

The only plausible explanation of the changes seen in 
Fig.[T] is the motion of KIC 11754974 in a binary system. 
The frequency changes can be interpreted in terms of the 
Doppler effect and transformed to radial velocities via the 
equation 



Vr, 



/rof 



where c represents the speed of light. These radial velocities 
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can be used to derive parameters of the spectroscopic orbit 
of the pulsating primary. It can be seen from Fig.[T] that 
the peak-to-peak radial velocity of KIC 11754974 amounts 
to about 16kms~^. 

In order to derive parameters of the spectroscopic or- 
bit, however, we proceeded in another way. Instead of using 
radial velocities calculated from frequencies, we decided to 
use the O — C diagram for the times of maximum light. 
The advantage of using the O — C diagram is that in this 
diagram the changes of period accumulate. For large orbital- 
to-pulsation period ratios, spectroscopic elements can be de- 
rived from the O ~ C diagram more accurately than when 
using radial velocities calculated from the Doppler effect. 

Since the amplitude of the strongest mode is several 
times larger than the amplitudes of the next three modes, 
we decided to use only the times of maximum light of the 
strongest mode. First, all terms with amplitudes higher 
than 40ppm except for /i and its harmonics were removed 
from the data. The residuals were divided into 5-20 d non- 
overlapping subsets. For each subset, a sinusoid with fre- 
quency /i and its harmonics were fitted in order to obtain 
amplitude and phase information. The time of maximum 
light was calculated from the phase of /i with the epoch 
adopted to be the nearest to the average time of observa- 
tion in a given subset. The resulting observed (O) times of 
maximum light were compared to the calculated (C) times 
of maximum light according to the following ephemeris: 

C{E) = 2454950.021179 + 0'?0611817444 x E, 

where E is the number of cycles elapsed from the initial 
epoch, in BJD using Coordina ted Universal Time (UTC; 
lEastman. Siverd fc Gaudill2010f ). As such, our Kepler times 
do not suffer the timing error discovered in 2012 Novembeijj. 

The Wide Angle Search for Planets (WASP) is a multi- 
site, multi-camera imaging system obtaining photometric 
data with an a ccuracy better than 1 per cent for objects with 
7 ^ V ^ 11.5 (jPoUacco et aLlbOOd ). Lower-quality data are 
obtained for stars with 11.5 ^ V ^ 15. Although the Kepler 
data have much higher precision and duty cycle, the WASP 
observations have time-spans extending back to years before 
Kepler came online. We utilised this longer time-base of the 
WASP data to confirm our O — C results, the values of which 
are plotted in Fig.[2l 

Using the more precise Kepler times of maximum light 
(for both LC and SC data) , we fitted the O — C parameters 
assuming an eccentric orbit. The parameters of the fit are 
given in Table[T] The meaning of the parameters is standard: 
To denotes the epoch of periastron passage; e, eccentricity; 
Lu, longitude of periastron; Ki, half-range of the primary's 
radial velocity variation; Porb, the orbital period; Po, the 
pulsation period; ai, the semi- major axis of primary's ab- 
solute orbit; and f{M), the mass function. The residuals 
from the fit are equal to 3.28 s for the LC and 1.47 s for the 
SC data. At least part of the residual scatter comes from 
the small-amplitude variation due to modes that were not 
subtracted from the data. 

The lower panel of Fig. [2] indicates the WASP data 
have large residuals (in time) when compared with 
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Table 1 . Parameters of the spectroscopic orbit derived from the 
O — C diagram. To is given in BJDutc- 



Parameter 



Value 



To 


2454999 ± 37 


e 


013 +0011 


^n 


102 ± 38 


K^ [km s-i] 


8.35 ± 0.04 


Porb [d] 


343.27 ± 0.34 


^o[d] 


0.0611817444 ± 0.0000000011 


a\ sini [AU] 


0.2634 ± 0.0013 


/(M) [Mo] 


0.02069 ± 0.00031 



the Kepler data. We tested whether the shift is in- 
trinsic by comparing two stars in which such shifts 
were deemed to be unlikely - we chose two relatively 
bright, short-period WUMa stars, namely: KIC 8554005 
and KIC 9392683 (1SWASPJ191805.45-H444115.4 and 
1SWASPJ190413. 72-^455657.6, respectively). We inves- 
tigated the strongest (second) harmonic in the Fourier 
transform of each star, and define a parameter DT = 
rmax(WASP) - Tmax (Kepler), where Tmax represents the 
time of maximum light closest to the mean epoch of a given 
dataset. For KIC 8554005 we determine DT = -5.8 ± 5.2 s 
and for KIC 9392683 DT = -13.7 ± 6.0 s, and note that 
DT is therefore zero to within 1-2 a. For KIC 11754974 we 
found DT = -34.2 ± 2.8s from aU data or -27.2 ± 2.1s 
if the first season is excluded. We therefore conclude that 
the shift for KIC 11754974 is intrinsic and probably occurs 
because the WASP and Kepler bands are different; the 
central wavelengths are ~550 nm for WASP and ~660 nm 
for Kepler. We note that WASP and Kepler use heliocen- 
tric and barycentric Julian date, respectively, and that 
contributes up to a couple of seconds to the offset. 

One can see that the orbit has a small eccentricity. 
The secondary is likely to be a low-mass main-sequence 
star. Assuming the primary's mass lies in the range 1.5- 
2.5 Mq, we get M2 = 0.42-0.58 M© for inclination i = 90° 
and M2 = 0.65-0.87 Mq for i = 45°. This indicates a K or 
early M-type secondary companion. Its contribution to the 
total fiux is therefore not larger than a few per cent unless 
the inclination is very low. This reduces the possibility of 
detection of variability of the secondary if present at all. 

Having calculated the orbital elements we have cor- 
rected the BJD times for short-cadence Q6 and Q7 data for 
the light-time effect in the binary system. Because in differ- 
ent quarters different detectors and apertures are used, it is 
expected that fluxes from different quarters are not directly 
comparable and need to be scaled. Using the amplitudes of 
the strongest modes we found that the amplitudes derived 
from Q7 data are larger by a factor of 1.0087 than those 
obtained from Q6 data. Therefore, the Q7 data were di- 
vided by this factor prior to combining. Subsequently, this 
high-duty cycle combined Q6+Q7 SC dataset was subject 
of time-series analysis (Section[4]), where we have used only 
these adjacent SC quarters to avoid aliasing difficulties. 



2.2 Frequency modulated stars 

The effect of the binary system on light arrival times 
is that of frequency modu lation for the pulsation modes. 
IShibahashi fc Kurt a (|2012l ) showed how a frequency multi- 
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Figure 2. The O — C diagram for tlie times of maximum liglit of the strongest mode of KIC 11754974 (upper panel). Filled and open 
circles denote Kepler LC and SC data, respectively. The open squares are the 2007-2010 WASP data. The lower panel shows the residuals 
from the eccentric orbit fit, shown as a continuous line in the upper panel. The parameters of the fit are given in Table[T] 



plet is induced around each pulsation mode as a result of the 
binarity. The frequency spacing can be used to infer Porb, 
and the amplitudes and phases of the frequency multiplet 
give To, e and f(M), provided the dataset is longer than 
Porb. .f{M) is calculable using the ratio of the amplitudes of 
the first orbital sidclobcs to those of the central pulsation 
frequency fShibahashi & Kurtz 2012). We used this method 
to check our binary parameters, and we will show that the 
results we obtain using this method are entirely consistent 
with those derived from the O — C diagram. We applied 
this methodology to the Kepler LC QO-13 dataset, span- 
ning 1153d (3.36 orbits). 

We started with the highest-amplitude peak in the 
QO-13 dataset, which lies at 16.344745 d~^ with an ampli- 
tude of 51.79 mmag. We extracted this peak and its orbital 
sidelobes and applied non-linear least-squares fitting rou- 
tines to improve the frequencies, amplitudes and phases. 
From the mean spacing of the multiplet we derive Po^b = 
343.66 ± 0.48 d, which is consistent with that derived via 
the O — C diagram. Given the slightly lower uncertainty on 
Porb from the latter method, we used the (O — C)-calculated 
value of Porb to force-fit the sidelobes to be exactly equally 
spaced. We chose the zero point in time, such that the side- 
lobes have the same phase with respect to each other, but 
are 7r/2 out of phase with the central frequency, thus demon- 
strating the triplet is orbital in nature. Table[2] displays the 
results for this frequency and three others, each of which 
shows a frequency triplet split by the orbital frequency. 

The mass function is calculable as: 



/(mi, 7712, 



1+1 



+ A^ 



Ao 



^o'rb 2vrG 



where the amplitudes of the central, higher and lower fre- 
quency components of the triplet are denoted Ao, A+i and 
A-i, respectively; Pose is the oscillation period; and all quan- 
tities are in SI units. Via this method we find f{M) = 



0.02069 ±0.00036 Mq, which is exactly the value calculated 
via the O - C method. Eq. 23 of lShibahashi fc Kurt j (|2012l ) 
allows one to calculate ai sini, too, as 



ai sm I 



27r 



Application of this formula yields aisini = 0.2629 ± 
0.0051 AU, in good agreement with the O — C value. 

That such information is discernible from the Fourier 
transform highlights the importance of having continuous 
long-term observations available for 5 Set stars, whose high- 
amplitude, high-frequency pulsations produce sidelobes with 
high signal-to-noise (S/N). 

We also looked for eclipses in the SC data at the times 
when the stars are aligned to the line of sight. The possibil- 
ity of finding such a shallow eclipse in the residuals of the 
light curve of a binary with this long a period, even when 
our multi-frequency fit was subtracted, was very small. Al- 
though we did not find an eclipse, for such a wide system this 
provides only a slight constraint on the orbital inclination 
(an upper limit that is near 90°). 



3 SPECTROSCOPIC OBSERVATIONS 

We obtained three spectra of KIC 11754974 with the FIES 
spectrograph at the 2.5-m Nordic Optical telescope (NOT). 
The star is too faint for high resolution observations with 
these facilities so to obtain spectra of a higher S/N we ob- 
served with medium resolution (R=25000). Even still, ob- 
serving conditions rendered the S/N ratio of the obtained 
spectra low - the S/N ratio of the averaged spectrum, esti- 
mated from the noise level of reduced ID spectra, is about 
45 in the spectral region of Ha and about 40 near H/3. The 
difficulties in determining the right continuum level limit the 
accuracy of our determinations, which are provided to offer 



© 0000 RAS, MNRAS 000, 000-000 



KIC 11754974 5 



Table 2. A least-squares fit of the frequency triplets for tiie four liigiiest-amplitude modes. Tlie sidelobes are fixed at frequencies of 
Voac ± t'orb- Tile zero point in time, to = BJDutc 24555331.39590, was chosen so that the sidelobes have equal phases and the difference 
in phase between the sidelobes and central peak is 7r/2. Column 4 thus shows that the difference in phases of the sidelobes with each 
other is zero within the errors, and Column 5 demonstrates their phase difference from the central peak is 7r/2 within the errors. Column 
6 contains the amplitude ratios of the sidelobes to the central frequency. Since these values are small compared to unity, they are 
approximately equal to a, the amplitude of the phase modulation. Column 7 shows the expected theoretical outcome that the ratio of a 
to the oscillation frequency is the same for all triplets. We note that for the first triplet, the amplitudes of the sidelobes arc only equal 
at the Act level, and that the errors provided arc overestimated because there are still many frequencies left in the data. 

a = (A+i + A_i)/Ao a/i^osc 

xlO-^d 



frequency 


amplitude 


phase 


4>+i — (t>-i 


(0+1 -'t>-l) - 


d-i 


mmag 


radians 


radians 


radians 


16.341823 


3.964 ± 0.048 


-1.6554 ±0.0120 






16.344745 


51.792 ±0.049 


-0.0545 ± 0.0009 


-0.011 ±0.017 


-1.606 ± 0.008 


16.347667 


4.148 ±0.048 


-1.6663 ±0.0115 






21.396067 


1.116 ±0.048 


0.9328 ± 0.0427 






21.398989 


11.046 ±0.049 


2.5642 ± 0.0044 


0.012 ±0.060 


-1.626 ±0.030 


21.401911 


1.140 ±0.048 


0.9443 ± 0.0418 






20.904493 


0.766 ± 0.048 


0.0139 ± 0.0623 






20.907414 


8.320 ± 0.049 


1.6598 ± 0.0059 


0.045 ± 0.086 


-1.623 ± 0.043 


20.910336 


0.803 ± 0.048 


0.0589 ± 0.0595 






20.940632 


0.584 ± 0.048 


2.1123 ±0.0818 






20.943554 


5.890 ± 0.049 


-2.5581 ± 0.0083 


-0.003 ±0.114 


-1.614 ±0.058 


20.946475 


0.603 ± 0.048 


2.1090 ± 0.0792 







0.157 ±0.003 



9.58 ±0.16 



0.204 ±0.012 9.54 ±0.58 



0.189 ±0.016 9.02 ±0.78 



0.202 ±0.023 9.62 ±1.10 



some constraints on fundamental parameters, rather than 
as definite determinations. We stress that higher resolution, 
higher S/N spectra are still needed. 

The Balmer-line profiles are good indicators of effective 
temperature for A-type stars, thus we have compared ob- 
served and synthetic profiles of the Ha and H/3 lines. The 
synthetic calcula tions of Ba l mer p rofiles were done using the 
SYN TH code by IPiskunovl (119921 ) and model atmospheres 
from iHeiter et al. (120021') with sol ar metallicity (as deter- 
mined bv lAnders fc Grevessdl 19891 ). The source of the spec- 
tral line list was th e Vienna Atomic Line Database (VALD; 
iKupka et al.l 1 19991 ) . Later analysis of metal lines demon- 
strated metal deficiency with respect to solar abundances, so 
we switched to working with low metallicity models ([M/H] 
= —0.5). For Ha and H/3 the best fits for observed and syn- 
thetic profiles were obtained for this metallicity model, yield- 
ing Tcff = 7000 K and 6800 K, respectively. Considering pos- 
sible errors in continuum determination of several percent, 
we estimate from the Balmer lines that T^s = 7000 ± 200 K. 

With our spectral material we cannot make reliable 
determinations of other parameters, namely logg and the 
microturbulence velocity Wmic- We calculated a number 
of synthetic spectra keeping Toff and one other parame- 
ter fixed, (thus permitting the third variable to change) 
to find the best fit of the metal lines. After several it- 
erations we accepted a model atmosphere with Tcff = 
7000 K, \ogg — 3.6 (cgs) and did calculations for micro- 
turbulence velocity resulting in Wmic = 2kms~^. This logg 
value is only loosely constrained and we considered that 
the star shows HADS-star characteristics and is proba- 
bly more evolved. However, the logg value obtained is 
lower than we expected, gi ven that cool S Set stars tend 
to pulsate in low overtones (|Breger fc Bregmanlll975l ). and 
a low overtone mode at 16 d~^ (the frequency of our 
highest-amplitude mode) implies an evolutionary status 
much closer to the ZAMS. We compared the spectrum with 
the spectrum of S Set itself, taken from the Elodie archive 
(|http : //atlas . obs-hp . f r/elodie/ [) . The shapes of the H/3 



profiles are very similar in these stars. However, the spec- 
trum of S Set is richer in metal lines. 

The weak spectral lines in KIC 11754974 are lost in the 
noise. In the averaged spectrum we identified lines of Fel, 
Fell, Till, Mgl, Mgll and some others. For identification 
we calculated a synthetic spectrum using our best fitting pa- 
rameters. From several lines including Fell A5018.440, Mgl 
A5172.684, Mgl A5183.604 we determined a projected ro- 
tational velocity: vsini = 31 ± 2kms~^. A more cautious 
interpretation, taking into account more lines (with a poorer 
overall fit) is vsini = 25 ± 7kms~^ - given the low S/N of 
our spectrum, we adopt the more cautious vsini value. 

Using several lines we determined the abundances of 
iron: log N/Ntot = -5.8 ± 0.3 (while solar is -4.54Jf|, and 
for magnesium: log N/Ntot = —5.9 ± 0.3 (solar is —4.46). 
We note that fitting Ha and H/3 with [M/H] = -1.0 yielded 
negligible differences from the [M/H] = —0.5 models. Many 
other chemical elements in the spectra also show significant 
deficiencies but the lines are blended or below the detection 
threshold and cannot be tested. There is a problem in that 
many lines of Fe I that are not visible in the observed spectra 
appear in synthetic ones and require yet lower abundances 
for best fitting than we obtained from available lines, the lat- 
ter mostly belonging to Fe II. One must decrease the effective 
temperature to reduce such a disequilibrium. At least one 
strong available line of Fel 4404.750 A shows similar abun- 
dances to the Fe II lines. This problem needs higher quality 
spectra for further analysis. Stromgren photometry is de- 
sirable for independent determinations of T^a and log g pa- 
rameters. Our parameter determination from spectroscopy 
is given in Table[3l 

With these restrictions in mind, verification of the low 
metallicity was provided by a semi-automated routine for 
476 lines from the VALD database in the 5000-5200 A wave- 
length range. The best fit [M/H] value using x^ minimisation 



as indicated in logarithmic form in the VALD. 
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Table 3. Atmospheric parameters determined from spectroscopy. 



Tcff 
K 


log 9 

(cgs) 


[Fe/H] 
dox 


Dsin j 
kms~^ 


kms"-*- 


7000 ± 200 


3.6 ±0.3 


^s+o-^ 


25 ±7 


2 



was — 0.58dex, when the logg value was held at the deter- 
mined value of 3.6. The [M/H] value is not the [Fe/H] value, 
but is close to it. 

An independent check on our spectroscopic TcS can 
be obtained from the InfraRed Flux Method (IRFM) 
(jBlackwell fc Shallidl 19771 ) . There is no sign of any interstel- 
lar NaD lines in the spectrum, so reddening is expected to 
be negligible. Broad-ba nd APASS photometry from UCAC4 
IjZacharias et al.ll2012l ) and 2MASS were used to estimate 
the observed bolometric flux. The IRFM was then used to 
determine Tcff = 7110 ± 150 K, which is consistent with the 
Tcff from the spectral analysis. 

As another check, we also calculated a quartic poly- 
nomial fit in (<;' — r') colour to calculate a temperature of 
7170 K, but the uncertainty on th is value is large at about 
±250 K jUvtterhoeven et al.llioTTI '). 

We used the spectrum to discern the star's kinemat- 
ics, finding a velocity shift of — 300kms~^, in agreement 
with lBalona fc Nemej (|2012l '). We expect that the 16-km s"^ 
radial velocity variations inferred from the light curve (cf. 
Fig.[T]) would be detectable from spectroscopic observations 
that were suitably spaced in time. 

Underabundances of Fe-peak elements by factors of 
10-100 are characteristic of the A Boo stars, which com- 

J irise only ~2 per cen t of objects in the 5 Set spectral region 
Paunzen et al.ll2002l ). If this star is determined through a 
full abundance analysis to be of the A Boo class, it will be the 
most intensely photometrically studied A Boo star to date 
- only eight spectroscopic binary systems with suspected 
A Boo components are currently known l| Paunzen et al.l 
[2013). With our low S/N spectrum it is not possible to de- 
termine whether C, N and O are normal, as in the A Boo 
stars, but many elements seem underabundant with respect 
to the Sun. The observed abundances are entirely consistent 
with the classification o f KIC 11754974 as an SX Phe star by 
iBalona fc Nemed l|2012l ). 



4 PULSATION CHARACTERISATION 

4.1 The data and preliminaries 

In this paper we analysed the short cadence Q6+Q7 light 
curves from KASOQj (the Kepler Asteroseismic Science Op- 
erations Center) for the HADS star KIC 11754974, to which 
we make corrections to the times to account for the light- 
time effects associated with being in a binary system as men- 
tioned in §[21 We used the Pre-search Data Conditioned data 
of the least-squares pipeline (PDC-LS), but made checks 
with the Simple Aperture Photometry (SAP) data on which 
only basic calibration is performed in the data processing 
pipeline. 



The PDC data are useful for their satisfactory re- 
moval of instrumental signals, including, but not limited to: 
heating/cooling surrounding safe-mode events, cosmic ray 
events, and return to scientific operating focus and temper- 
ature after monthly downlinks (Kepler Data Characteris- 
tics HandboolQ). The PDC data do not necessarily treat 
all jumps and outliers. For KIC 11754974 a manual check 
prior to analysis showed no clear outliers in need of removal; 
after fitting 166 frequencies, 103 'outliers' of the 256 514 
data points were manually discarded, 100 of which were at 
least 1 mmag from the fit, and the independent frequencies 
were thereafter improved with non-linear least-squares fit- 
ting. While in a limited number of cases the PDC data mod- 
ify some stellar variability, known cases are documented in 
the data release notes available on the Kepler Guest Ob- 
server websitqj, and the benefits of using PDC data over 
SAP data outweigh the detriments in this case. Further dis- 
cussion on the implications to asteroseismology of choos- 
ing which of th e SAP or PDC data to use can be found in 
lMurphvll|2012al lbl). 

We also note that the variability of KIC 11754974 was 
already discovered fr om the All Sky A utomated Survey 
(ASAS) observations (|Pigulski et al.ll2009l . star #143 in the 
catalogue). Due to the quality of the data, only the main 
mode was detected; we did not use those data in this inves- 
tigation. 

The data were imported into PERIOD04, for its useful 
graphical user interface, convenient least-squares and dis- 
crete Fourier transform functions, and tools for frequency 
extraction with prewhitening. A review of the functions of 
PERIOD04 has been published by the program's creators 
(jLenz fc Brcger 2004 ) , and a more comprehe nsive view of its 
capab ilities can be found in the user guide (jLenz fc Bregeij 
l2005h . 

The corrected, edited Q6.1-Q7.3 dataset con- 
tains 256 411 points taken from (BJD) 2455372.439 to 
2455552.558 (180.12 d), having a 97 per cent duty cycle. 



4.2 The light curve 

We extracted 166 frequencies that are an excellent fit to the 
light curve. We recalculated our fit both after outlier ex- 
traction, and after removing a cubic spline fit. The latter 
was done by calculating the average residuals at 2-d inter- 
vals, interpolating with a cubic spline, and subtracting the 
smoothed curve from the original data. Semi-regular cycles 
(Fig.El lower panel) are evident on different timescales as a 
result of beating between frequencies. 

The peaks in brightness occur over a greater inten- 
sity range than the troughs (Fig.[31 upper panel). This 
kind of light curve shape is typical for HADS stars and 
can be thought of a s saturation of the driving mechanism 
(|Balona et al.l I2OIII ). It has also been observed in many 
7 Dor stars, and is common in RRab (RR Lyrae) vari- 
ables. In particular, it is common in white dwarfs that show 



many combination frequencies in their pulsation spectra. |W 
((2OOII) described how in white dwarfs, which pulsate in g- 
modes, the convection turnover time is inferred from the 



|http : //kasoc ■ phy s ■ au ■ dk | 



http : //archive ■ stsci ■ edu/kepler/manuals/Data_Characteristics .pdf | 



http : //keplergo ■ arc ■ nasa ■ gov/Documentat Ion ■ shtml | 
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Figure 3. 1.1-d segment of the light curve of KIC 11754974, showing obvious beating. Flux maximum occurs over a greater range in 
brightness than does flux minimum, as indicated by the fact that <5max > 'Jmin- The lower panel indicates this is not just localised to a 
short time segment. 



phase of a combination frequency relative to its parents, and 
how the amplitude difference between the sums and differ- 
ences in combination terms (e.g. /i +/2 and /i — /2) can give 
the thermal constant of the stellar convection zone at equi- 
librium. However, in a n investigation into combination fre- 
quencies in S Set stars, iBalonal (|2012h concluded that there 
is no information in the relative amplitude of a combination 
frequency that might be useful for mode identification in 
these p-mode pulsators. 



4.3 The region of independent frequencies 

In KIC 11754974 all of the high-amplitude independent fre- 
quencies are found in the range 16 to 25 d~^ (a full list of 
extracted frequencies with their corresponding amplitudes, 
along with identifications, is given in Tabled). Outside of 
this range, most frequencies detected can be explained as 
combination of dominant frequencies in the independently 
excited region. Such is also the case for the star 44Tau, 
on which the literature is extensive. In KIC 11754974, in- 
dependent peaks outside the 16 to 25 d~^ range have low 
amplitudes, typically of order 200 /imag. 

In HADS stars the period ratio of the two most dom- 
inant modes offers easy m ode identification when that ra- 
tio is 0.77 IjPorettil 120031 . for example), as this is the 
theoretical period ratio of the first radial overtone mode 
(hereafter Fl) to the fundamental radial mode (FO). The 



period ratio can be plotted against the base- 10 loga- 
rithm of the period of FO, logP(FO), in a so-called Pe- 
tersen diagrarrlj, offering diagnostic inform ation on the star. 
[Petersen &c Christensen-DalsgaardI l|l996r) discuss the dia- 
grams in detail for double-mode HADS stars of different 
metallicities in their figure 3. They showed that lower values 
of the P(F1)/P(F0) ratio are found for metal rich stars and 
stars closer to the ZAMS, but spectroscopy of KIC 11754974 
indicates that this star has neither property. Hence the pe- 
riod ratio of its two highest-amplitude modes, /i and /2, 
having P(/2)/P(/i) = 0.7638, is much below the anticipated 
value (> 0.770) from the metallicity-calibrated Petersen di- 
agrams. We can thus argue that /i and /2 are not the FO 
and Fl modes. We also note that the high-amplitude inde- 
pendent frequencies around 21 d^^ (Figs.U&O, particularly 
/a and /4, could each potentially be the first radial overtone 
based on period ratios, though not all can be ^ = modes, 
and none exactly gives the anticipated ratio. 

The high amplitude of /i alone does not necessar- 
ily imply it is the fundamental radial mode, either; in 
FGVir whose pulsation amplitudes are comparable to 



^ The diagrams were first discussed by iPetersenI 1 119781) , but it 
was Art N. Cox who brought the term 'Petersen diagram' into 
general use through his conference talks and eventually in his 
journal articles. 
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this star, the fundamental radial mode was not found 
to b e the highest in amplitude of the observed frequen - 
cies iJGuzik. Bradley fc Tem plcton'200^:'Breger ct al.'2005[). 
4CVn is another exam ple ( Breger, Davis fc Dukes 20081 : 
ICastanheira et al.l2008r ). The high amplitude of /i does have 
other consequences: its window pattern dominates the entire 
analysed spectrum, even far from its frequency of 16.345 d~^. 
Peaks arising solely from the spectral window of /i have am- 
plitudes of 0.9 and 0.5 mmag at 50 and 100 d~^ respectively. 
To make other peaks with lower amplitudes visible, Fig.|4] 
shows a schematic amplitude spectrum with all window pat- 
terns removed. The five highest amplitude peaks in the range 
16 to 25 d~^ are all independent modes (note that the third 
and fourth highest peaks in that range are indistinguishable 
in Fig.|3|. Other peaks with high amplitudes found outside 
this region are combinations of these frequencies. 

Only once the five main independent modes are removed 
is it possible to see the high frequency density in the Fourier 
transform. Fig. fallows a clearer view of the lower-amplitude 
peaks surrounding the main pulsation frequencies. In this 
frequency-range independent frequencies are common, and 
few frequencies appear to be combinations. The spacing of 
0.036 d~^ between /a and /j is found repeated in combina- 
tion frequencies all over the spectrum, with different coeffi- 
cients of /i and /2 added to them. 



4.4 Quintuplet 

The Fourier transform of the light curve of the star features 
a quintuplet, with a central frequency at 20.242 d~^, and 
two lower-amplitude companions on each side with a mean 
separation of 0.218 d~^, shown in Fig. |6] (upper panel). We 
assume that the quintuplet is rotationally split. The exact 
morphology of the quintuplet depends on excitation and ori- 
entation, but unless one can safely assume components are 
excited to the same amplitude, the amplitude ratios can- 
not be used to constrain the star's inclination. The separa- 
tions are not perfectly identical between components, nor 
do we expect them to be for rotationally split quintuplets - 
exact frequency spacing is not e xpected unless there is fre- 
quen cy locking due to resonance (jBuchler. Goupil fc HansenI 
1 19971 ). The need for high-dispersion spectroscopy or multi- 
colour photometry is highlighted in this mode identification 
attempt - better spectroscopy can more tightly constrain 
both V sin i and the fundamental parameters required to im- 
prove the models (described in Section[Sl), allowing more cer- 
tain mode identification. 

In other case studies, peaks due to rotation have been 
found at low fre quency with highly-significant amphtude s 
(e.g. KIC 9700322 iBreger et al.ll201]] : lGuzik fc Bregeill201lh . 
Hence we searched the Q6 and Q7 SC data for a peak at low 
frequency that might correspond to the rotation frequency. 
There is no evidence for a peak at the anticip ated value, nor 
at half or double that value (cf. lBalonall2011h . In fact, there 
is no evidence for a peak due to rotation below ~0.8d~^, 
which has implications for our modelling discussion ij |5.3l 
To check that such a peak was not somehow removed by 
either our spline fit or the PDC pipeline, the SAP data were 
also searched for the same peak, but no peak was found. 

Another quintuplet occurs at ~ 36.59 d~^ (Fig.fBl lower 
panel), but this is a manifestation of the main quintuplet 
in combination with /i . This is confirmed by inspecting the 
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Figure 7. Schematic diagram of the many quintuplets in the 
data. The main quintuplet at 20.242 d~^ forms other quintuplets 
throughout the Fourier spectrum in combination with /i . One can 
see the similar amplitude ratios of the quintuplets as the shape 
of the quintuplet is maintained in each combination. 



amplitudes of this second multiplet in light of the first - not 
only are they smaller, but they retain similar amplitude ra- 
tios between peaks within the multiplet. Similar quintuplets 
can be identified at 3.90 d~^, arising from /quintuplet ^ fi, 
and at 12.45 d~^ from 2/i — /quintuplet. The frequency dis- 
tribution and the amplitudes of the quintuplets are shown 
in Fig.[7] 



4.5 Combination frequencies 

The frequencies in Table|4]have been extracted down to an 
amplitude level of 100 /^mag for the region 0-70 d~^, and 
to 20/imag for 70-100 d~^ for closer examination of pat- 
terns. The region 70-100 d~^ shows the amplitude level at 
which pulsation-mode candidates can be identified. Kepler 
SC data allow high frequency analysis, and the high preci- 
sion of the satellite's photometry, along with its position out- 
side of Earth's atmosphere, afford low noise levels (< 5 /imag 
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for the region 70-100 d~^). The combination of the two make 
the detection of high order frequency combinations possible. 

The tolerance criterion for each calculated combina- 
tion frequency has been chosen as 0.001 d~^; to be iden- 
tified as a combination, the calculated frequency must agree 
with the observed frequency within this value. This toler- 
ance is strict by some definitions - for the dataset analysed, 
1/r = 0.0056 d~^ and the full width of a peak at half max- 
imum (FWHM) is 0.0067 d. A stricter tolerance reduces the 
probability of false positives for frequency combinations, and 
has been made necessary because of the high number of fre- 
quencies extracted. 

In addition to meeting the tolerance criterion, the com- 
bination must also be physically sensible - high coefHcients 
and many combining independent frequencies require in- 
creasingly unlikely interactions between modes. Hence both 
the coefficients of the independent frequencies and the num- 
ber of independent frequencies must be small. Sensible com- 
binations typically have neither of those numbers greater 
than three, but no hard-and-fast rule was applied. In the 
region above 70 d~^, the coefficients have to be quite high, 
especially since /i is only 16.345 d~^. We note that some 
of the remaining 'independent' frequencies without identifi- 
cations in TableUcan be explained with physically sensible 
combination frequencies if the tolerance is relaxed. 

PERIOD04 is capable of testing for frequency combina- 
tions to a given toleran c e crite rion; such a method was em- 
ployed by iBreger et al.l (|2005l ) in their analysis of the star 
FGVir, where their tolerance criterion was 0.001 d~^. We 
used the same tolerance criterion and ran simulations to 
determine the number of accidentally-matched frequencies 
one can expect. Using the frequencies /i to /s, allowing co- 
efficients between —1 and 6, and allowing all five frequen- 
cies to be involved in the combinations simultaneously (i.e. 
there is much more freedom than in the combinations in 
Table[3)), we find only 1-2 frequencies in the range 70 to 
100 d~^ to be accidentally matched. Hence almost all com- 
binations presented in Tabled are probably genuine. The 
two frequencies with the greatest uncertainties lie at 91.377 
and 95.939 d~^, with errors of 0.001 d~^, which is the same 
variation as the combination frequency tolerance. These are 
the lowest-amplitude frequencies extracted, and are the two 
prime candidates for mis-matched frequencies. 

Combination frequencies are even detected at beyond 
100 d"""^ (not shown). Particular caution has been made in 
identifying high frequencies as combinations, especially be- 
cause they often contain more than two independent fre- 
quencies and coefficients of 3, 4 or even higher. There is 
no doubt however, that these are indeed combinations. Sev- 
eral patterns of frequencies repeat themselves, separated by 
/i (Fig. [S]), in a way that could only happen dependently. 
The frequencies of these peaks are entirely calculable given 
the frequencies of the main independent modes, and at yet 
higher frequencies they recur with diminishing amplitudes 
until ~140d~^ at which point their amplitudes become in- 
distinguishable from noise and the patterns become less dis- 
cernible when peaks are missing. 

4.6 Validity of extracted frequencies 

One could continue to extract frequencies down to lower and 
lower amplitudes, but this begs the question of where to 
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Figure 8. Repetitive patterns are caused by adding /i to other 
combinations of frequencies, hence panels A&D, B&E, and C&F 
are each separated in frequency by /i = 16.345 d~^. The fre- 
quency range displayed by each is shown in the top panel. The 
first peak in panel D is 4/2; 4/2 — /i is too low in amplitude to 
be distinguishable from noise, hence it is not seen in panel A. 
The cause of relative amplitude changes of neighbouring peaks 
following the addition of /i is unclear, but the absolute changes 
in amplitudes in adjacent panels (note the change in axis scale - 
units are /xmag) is a result of having a greater number of combi- 
nation terms. 



stop. lOO/xmag was chosen as a sensible limit, to extract the 
important frequencies for mode identification, and also for 
clarity: mode splittings become harder and harder to spot 
when the density of peaks increases, as happens when one 
lowers the extraction limit. On occasion, combination fre- 
quencies with amplitudes below 100 /^mag were extracted, 
such as when they are expected as part of a multiplet. 
There are additional complications to extracting frequen- 
cies to lower amplitudes. With a 180-d dataset, frequency 
resolution is 0.006 d~^. One cannot extract frequencies more 
closely spaced than this. Moreover, nothing is gained by con- 
tinuing extraction to lower amplitudes - combination fre- 
quencies become less reliable as the risk of false positives 
rises sharply with the number of extracted frequencies, mode 
identification becomes more difficult when mode density in- 
creases, and the low amplitude frequencies become 'noise' in- 
terspersed with the meaningful high-amplitude peaks. Fig. [9] 
shows the residuals left after pre- whitening all frequencies in 

Tabled 

The Kepler Input Catalogue (KIC; iBrown et aLll201ll ) 
provides information on the amount of light contamina- 
tion appearing in the data. The value of contamination for 
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Figure 9. Residuals left behind after extraction of 166 frequen- 
cies, to amplitude limits 100 /^mag for ^ / ^ 70 d^^, and 
20Mmag for 70 < / s£ 100 d"!. 



KIC 11754974 is 0.009, hence less than 1 per cent of light 
incident on the CCD pixels attributed to this star can be 
coming from background stars. As far as we can tell from 
looking at Kepler Full Frame Images, 2MASS and DSS im- 
ages, KIC 11754974 is safe from contamination. This makes 
negligible the possibility that a background star has a high- 
amplitude pulsation mode that is picked up and interpreted 
as a pulsation mode for KIC 11754974. We know the star is 
in a binary system, but with a star significantly less massive 
than the Sun, contributing little hght, and not expected to 
generate frequencies in the periodogram at the amplitudes 
and frequencies observed. 

It may be possible to check whether a frequency is 
a part of a combination or is an independent frequency 
by examining the relationship between phase and fre- 
quency for the combination terms. In a recent paper on 
KIC 9700 322, in which many combination frequencies were 
detected, iBreger et al.l (l201l|) found that the higher fre- 
quency modes have greater relat ive phases. The theory 
is bas ed on that developed by iBuchler. Goupil &: HansenI 
(Il997|) who described nonlinear resonances mathematically, 
and lBalonal l|2012h investigated further. A combination term 
with frequency / = ni/i +nj/j should have a relative phase, 
4>r = 4>c — {ni(j)i + nj(/!)j), where (j>c is the phase calculated 
by least-squares fitting. We find that (pr decreases with fre- 
quency once phase folding about 27r is taken into account - 
we shifted groups of points by multiples of 27r in phase until 
they lay roughly in a straight line, then made 27r corrections 
for individual points lying more than n away from the fit in 
an iterative way. All combination frequencies listed in Ta- 
bleware plotted in Fig. 1101 Those 'combination frequencies' 
that do not fit the correlation between (jtr and frequency 
in the figure might actually be independent frequencies, in- 
stead. On the other hand perhaps the relationship is a purely 
mathematical construct. The points lying furthest from the 
fit are 3/2 and 4/2, which given the nature of the other 
combination frequencies in Table[4] would be expected to be 
present and with high S/N. 

The correlation has been discussed by iBalonal l|2012h , 
who showed that for two interacting modes the phases of 
the harmonics of a given mode decrease linearly with in- 
creasing harmonic order (i.e. frequency). The observational 
evidence was not discussed in light of the precision with 
which frequencies, amplitudes and phases can be determined 
with continuous space-based observations. The effect of such 
small error bars is to reduce the significance of the observa- 
tional evidence substantially - the point with the median de- 







. . , 












' 








2 





-^■-^-^ 


fc. 
















- 


OJ 






« 


V *• 1. , 














~. 




-5 


r 






'i- • 












- 






































». 






































-lU 


- 












,1 


*■ 


>. 


•- 


Pi 


















• 




-., 



20 



40 60 

Frequency {d" ) 



100 



Figure 10. Relative pliases are negatively correlated with fre- 
quency for the combination frequencies. Those points located far 
from the others might be interpreted as independent frequencies 
instead of combinations. Error bars for phases are shown, but are 
generally much smaller than the plot symbols. Phases and errors 
are the formal least-squares values and were calculated with re- 
spect to the mid-point of the combined Q6 and Q7 SC dataset: 
BJDuTC 2455462.49792. 



viation lies over 27 ct away from the fit. The median was cho- 
sen because taking the mean is biased towards those points 
with minuscule relative phase errors, whose distance from 
the fit can be many hundred a. The combinations giving 
rise to such small errors typically involve /i and /2, the 
calculated phase errors of which are as low as 10"'' rad. 
This further supports a rejection of the hypothesis of any 
meaningful correlation, because combinations like /i -I- /2 
are superiorly statistically significant in terms of S/N, and 
should consequently lend greater weight in this correlation. 
The reduced chi-squared parameter, Xrcdi is on the order of 
10^, clearly indicating an inappropriate fit. 

To establish whether the correlation need be rejected 
entirely we ran a calculation whereby the frequency combi- 
nations and their phase errors were kept fixed, but instead of 
calculating phases from least-squares we randomised them. 
We found that a correlation cannot be constructed by adding 
±n27r to the random phases, so while it is clear that in the 
relative phases and frequencies of KIC 11754974 some corre- 
lation does exist, we conclude the correlation is neither tight 
nor useful. 



5 ASTEROSEISMIC MODELS 

KIC 11754974 lies near the red edge of the 5 Set instabil- 
ity strip. At these cooler temperatures, convection plays 
an important role. Indeed, without Time-Dependent Con- 
vection (TDC) treatments, an actual red edge is not accu- 
rately predicted. Thus the introduction of TDC models was 
an important step forwa r d in the model l ing o f 5 Set stars 
(|Dupret et al.ll2004l2005l ). lMurphv et al.l (|2012D gave a suc- 
cessful first application of these models to an individual 5 Set 
star observed by Kepler. 

Rotation plays a big role in the evolution and oscil- 
lations of ^Sct stars |Goupil et all l2005l : [Rodriguez et al.l 
l2006al lbl: IPoretti et al.ll201ll ). even for KIC 11754974, which 
rotates slowly compared to the distribution of 5 Set-star 
rotational velocities. Furthermore, asteroseismic analyses 
with rotating models can assist with mode identification 
through the effects of rot ation on the mode couplings and 
frequency relationships llSuarez. Garrido fc Goupill l2006l : 
ISuarez. Garrido fc Movall2007^ . 
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No combination of both TDC and rotating models ex- 
ists at present, so our modelling approach consists of the 
two methods applied separately. This section explains the 
details of the models used, and the methodology followed. 



5.1 Description of the models 

Our TDC models use the s ame tools and selection method as 
that explained in detail in iMurphv et al.l (|2012l ) ; the codes 
used ar e the Code Liegeois d'evolution stellaire (CLES; 
IScufl aire et al.l l2008h for structure models and the MAD 
code (IDupretl 2001f) supplemented with Gabriel's treatment 
of TDC (JGabriejIigQel : ICrigahcene et al.l l2005l ) for non- 
adiabatic non-radial oscillation calcula tions. Convection was 
treated using Mixing- Length- Theory (|Bohm- Vitensd 1 1958I . 
MLT), and TDC runs throughout the evolution calculations. 

In order to theoretically characterise 

KIC 11754974, pseudo-rotat i nq s ei smic models (see 
ISoufi. Goupil fc Dziembowskil Il998l : ISuarez et all 120021 ') 
of the star were built. These equilibrium models were com- 
puted with the evolutionary code CESAM (Morel 1997 ) using 
frozen convection (FC) and taking the strongest effects 
of rotation into account. The latter is done by including 
the spherically-averaged contribution of the centrifugal 
acceleration, which is included by means of an effective 
gravity g^g = g — Ac{r), where g is the local gravity, r is the 
radius, and Ac{r) — | rQ^{r) is the centrifugal acceleration 
of matter elements. The non-spherical components of 
the centrifugal acceleration are included in the adiabatic 
oscillation computations, but not in the equilibrium models 
(details in lSuarez. Goupil fc Morell2006l '). 

Since we have no information about the internal rota- 
tion profile of KIC 11754974, i.e., about how the angular 
momentum is distributed in its interior, we have adopted 
the hypothesis of differential rotation caused by local con- 
servation of the angular momentum (in the radiative zones) 
during the evolution of the star. Similar rotation profiles 
have been found when analysing the evolution of giant stars 
including rotationally-induced mixing of chemical elements 
and transport of angular momentum IjMaeder fc Mevnetl 
[2OOJ) . The physics of the equilibrium models has been cho- 
sen as adequate for i ntermediate-mass A -F stars, and was 
described in detail bv' Suarez et al.l ([2009). 

Adiabatic oscillations w ere computed using the adia- 
batic oscillation code filou (|Suarej|2002l : ISuarez fc Goupill 
I2OO8I ). This code provides theoretical adiabatic oscillations 
of a given equilibrium model corrected up to second-order 
for the effects of rotation. These include near-degeneracy 
effects (mode coupling effects), which occur when two or 
more frequencies are close to each other. In addition, the 
perturbative description adopted takes radial variation of 
the angular velocity (differential rotation) into account in 
the oscillation equations. The oscillation spectra were calcu- 
lated from frequencies around the fundamental radial mode 
(for each model) and the cut-off frequency. This allows the 
presence of low-order g and p modes (mixed modes), which 
are generally present in S Set stars. Since non-radial oscilla- 
tions are expected, modes were computed from degrees i = 
(radial) to f = 3. Although the calculations are adiabatic, 
the effects of rotation alter the frequencies by an amount 
greater than the use of non-adiabatic codes. 



5.2 Methodology 

The low Tcff of the target suggested TDC models would pro- 
vide a more accurate first approach to the global properties 
of the star. We selected the best-fitting model by perturbing 
the atmospheric parameters from §[3] (Table[3ll to obtain the 
best match with the observed oscillation frequencies, adjust- 
ing the TDC parameters (namely the convection efficiency) 
accordingly. The properties of the best TDC model are pre- 
sented in TableO 

It is important to stress that the TDC models do not 
include rotation. Our next step was therefore to calculate 
equivalent (in terms of global parameters) FC models, incor- 
porating rotation. We considered that the error committed 
this way (i.e. TDC first) is much lower than through direct 
modelling with FC and rotation. 

The modelling with rotating models involves building 
a grid, with masses ranging from 1.4 to 2.0 Mq in steps 
of O.lMfi^, a convective effic iency ranging from 0.5 (pre- 
scribed by ICasas et all |2006| for A-F stars) to 2.0 (which 
gave the best fit in TDC calculations) in steps of 0.1, and 
an overshoot parameter dov in the range of [0, 0.3] in steps 
of 0.05. The surface rotational velocity of the models, Vcq, 
ranged between 20 and 30 kms~^ offering enough dispersion 
to analyse different configurations of mass, metallicity and 
age. In principle this gives splittings close to those observed 
for the quintuplet, but this was not constrained to be the 
case. In addition, it indirectly guarantees a dispersion in the 
inclination angle probed within the model grid. Metallicity 
was limited to [Fe/H]= —0.5 ± 0.1 dex, though metallicities 
closer to the solar values were investigated and offered a 
poorer match with the observed frequencies. This implies 
the entire star is metal poor, and that we are not just ob- 
serving a surface abundance anomaly, thus arguing against 
a A Boo classification (cf. debates in th e literature about the 
interior metallicity of A Boo stars, e.g. iMova et al.ll20ld ). 

Model selection involves searching through the grid, 
minimising the mean-square error function x^- For each 
mode, 



1 ^ 



- J'th, 



(1) 



is calculated, where fohs,i and i/th,; represent the observed 
and theoretical frequencies respectively. The total number of 
observed frequencies is represented by A*'. Calculations are 
made considering the mass M, the rotational velocity Vcq, 
metallicity, omlt, dov, and age as free parameters; pulsation 
amplitude information is not directly used. Importantly, all 
the frequencies are fitted simultaneously; assumptions like 
the "strongest modes" being radial or quintuplet member- 
ship were not used as constraints for the models. 

We did look separately for diagnostic frequency ratios 
in the data. We took models of three difi'erent masses and 
searched for ratios of 0.772 ± 0.009 (cf. discussion in §|4]3} 
among the independent modes. Fig.[TT]displays the four pos- 
sible frequency ratios observed, none of which satisfies the 
Toft and (rather loose) log g uncertainties, or indicates a par- 
ticularly promising identification for the first radial over- 
tone. We stress this only reduces the likelihood that FO is 
present with Fl; either could still be present without the 
other. 
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Figure 11. Rotational Petersen Diagram covering models with 
1.45-1.60 Mq. Coloured in red are the models that satisfy the 
Tgff and logg uncertainties (7000 ± 200 K; 3.6 ± 0.5); we used 
a conservative logg uncertainty, because atmosphere models do 
not include rotation to calculate an effective gravity. Also plotted 
(green dots) are frequency ratios (e.g. /1/4 indicates /1//4) that 
fall within 0.772 ± 0.009 (see text). Numbers in brackets in the 
legend denote the value of Qmlt a-nd the overshooting parameter, 
respectively. 



5.3 Discussion of the models 

In general, we found that the models with the lowest x^ 
were below 1.60 Mq. Moreover, those that matched either 
FO or Fl to the observed frequencies had masses in the range 
1.50—1.53 Mq, which are very similar to the TDC-predicted 
masses. The best-fitting rotating model identifies /i as the 
fundamental radial mode, and the quintuplet as components 
of a mixed (n, £) — ( — 1, 3) mode. The parameters of this 
model are displayed in Table[Sl 

That the quintuplet is identified as a mixed mode has 
strong implications. It directly affects the Ledoux constant, 
Cl, depending on how much g-mode character such a mode 
has. For lower radial orders, Cl can exceed 0.1 for 1 = 2 
modes. Therefore, we had to search a range of low frequen- 
cies to account for this factor in f jAAl but still no peak due 
to rotation was found up to twice the expected value (i.e 
between 0.218 and 0.436 d~^). 

The asymmetry of the quintuplet gives a proportion- 
ality constant, Dl, for the second order effect (centrifugal 
force) of rotational splitting, which we can use as a consis- 
tency check. For moderately-rotating 5 Set stars, oscillations 
can be expressed in terms of the perturbation theory as 



■ t^o + (Cl ~ l)mn + Dl 






which describe s the effects of rotation up to second order 
(adapted from iDziembowski fc Goodelll992l ). for a rotation 
frequency, f2, and observed and intrinsic (unperturbed) pul- 
sation frequencies, oj and loq, respectively. Under the as- 
sumption that the quintuplet is actually an i" = 2 mode, one 
gets Dl=2A5 in the case of Cl=0, or L>l = 1.8 for Cl=0.12. 
These are abo ut the expected values for a n (n, i) = (2, 
2) mode. (See IDziembowski fc Goodd Il992l for the exam- 
ple values). The Q-value of an (n, £) — (2,2) mode at the 



observed quintuplet frequencies is consistent with the as- 
sumption that the lowest independent frequency is near the 
radial fundamental mode. 

Thus we have the following scenario: the best rotating 
model indicates that /i is the fundamental radial mode, and 
inferences on the Q-value of the quintuplet if it is an (n, £) 
— (2,2) mode, would agree with that identification for /i. 
However, the model predicts that the quintuplet is part of 
an ^ = 3 mode instead. This is only an apparent disagree- 
ment. First, that an (n, £) — (2, 2) mode can generate the 
observed frequency quintuplet does not make its identifica- 
tion as part of an £ = 3 septuplet impossible, as the expected 
rotational splitting asymmetry of such a mode is expected 
to be quite similar. Second, differential rotation in the ra- 
dial direction is another possible explanation. It has been 
proven that the oscillations of low-order p and g modes are 
highl y sensitive to variations of the ro tation profile near the 
core (jSuarez. Garrido fc Goupilll2006l ). Further attempts at 
direct mode identification, through multi-colour photome- 
try, and further constraints for the models from high-S/N, 
high-resolution spectroscopy are required for further inves- 
tigation. 



6 CONCLUSIONS 

KIC 11754974 is a cool, high-amplitude 5 Set star with a 
large number of combination frequencies that show a pre- 
dictable pattern at high frequency. We discovered from 
the pulsations that the star is in a 343-d binary system; 
our seismic models show the primary (pulsating) member 
is 1.53 Mq, making the secondary 0.63 Mq for a model- 
determined i = 47°. No direct contribution from the sec- 
ondary is seen in the spectrum or the light curve. 

Compared to most A-type stars, KIC 11754974 is a slow 
rotator, which is typical of HADS stars, and allowed us to 
construct perturbative rotating models. The models indicate 
the star is metal-poor, in that solar-metallicity models were 
a poor fit to the observed frequencies. This agrees with lit- 
erature classifications that this is a Population II (SXPhe) 
star, and argues against a surface abundance anomaly model 
like the Am stars or A Boo stars. 

The location of this star on the rotating Petersen dia- 
gram argues against more than one of the high-amplitude 
/1../5 modes being radial. We cannot confirm whether the 
quintuplet found in the data is rotationally split - direct 
mode identification is required for this. 

Combination frequencies are present all across the spec- 
trum, arising from five independent modes that are all con- 
fined to a small frequency region. Outside the frequency re- 
gion 16 < / < 25 d^^ independent modes with amplitudes 
above 500/imag are non-existent, and most frequencies are 
combinations of some form. 

The star has a lot of asteroseismic potential, being 
near the red-edge of the 5 Set instability strip where the 
pulsation-convection interaction is under study, being a Pop- 
ulation II star with precise and sensitive photometry, and 
being a high amplitude pulsator for which mode interaction 
is also particularly important. Here we have scratched the 
surface with our models, but there is great potential to learn 
much more. We intend to study the star further, with mul- 
ticolour photometry and higher-resolution spectroscopy to 
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identify the pulsation modes and constrain the atmospheric 
parameters. Despite being faint, the star is a promising tar- 
get for learning about mode selection in these pulsators. 
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Table 4. Multifrequency extraction and identification for KIC 11754972. The absence of combinations between 16 and 25 d~^ is evident; 
the quintuplet and /i to /s all lie in this region. Frequencies are only given an individual identification if they are of particular significance 
or form a combination. Errors on frequencies are generally smaller than 2xl0~*d~^, and errors on amplitudes are all under 4/imag. 
For /i to /5, from which combinations are calculated, the frequency uncertainty is 10~®d~^. * denotes quintuplet membership, f this 
frequency is not fully resolved from /2 in the SC data. 



frequency 


amp. 


ID 


frequency 


amp. 


ID 


frequency 


amp. 


ID 


d-l 


^mag 




d-i 


fmiag 




d-i 


^mag 




16.34474 


63329 


/i 














21.39898 


15940 


/2 














20.90740 


11653 


/3 














20.94355 


8229 


/4 














22.66020 


2839 


/5 














0.4554 


270 


= /2 - /4 


20.4552 


603 


= /ii* 


42.7980 


946 


= 2/2 


0.4916 


865 


= /2 - /3 


20.6552 


422 


= /l2* 


43.5676 


122 


= /3 + /5 


0.7723 


202 


= /7 - /l 


21.0537 


216 




44.0592 


389 


= /2 + /5 


1.0705 


315 


= /l - /6 


21.3973 


257 


t 


44.4716 


107 


= 4/1 - /3 


1.2612 


619 


= /5 - /2 


21.4282 


274 




47.8522 


103 


= 3/2 - /l 


1.3804 


178 


= /2 - /9 


21.7251 


458 


= /l4 


49.0342 


2677 


= 3/1 


1.7528 


250 


= /5 - /3 


22.5325 


270 




53.5969 


1230 


= 2/l+/3 


3.4387 


202 


= h-fi 


22.8592 


107 




53.6330 


828 


= 2/1 + h 


3.6738 


291 


= h-h 


22.8860 


205 




54.0885 


1466 


= 2/l + /2 


3.8976 


334 


= ho - h 


23.7744 


354 


= /l5 


55.3497 


235 


= 2/l+/5 


4.1040 


141 


= /13 - /l 


24.3270 


130 




58.1595 


246 


= /l + 2/3 


4.1105 


161 


= /ii - h 


25.4701 


257 


= 2/3 - /l 


58.1957 


174 


= /i + h + h 


4.3108 


76 


= /12 - /l 


25.5062 


229 


= /3 + /4 - /l 


58.6511 


292 


=h+h+h 


4.5627 


4625 


= h-h 


25.6208 


242 




58.6873 


191 


=h+h+h 


4.5988 


2234 


= fA-h 


25.9616 


403 


= /2 + /3 - /l 


59.1427 


705 


= 2/2+/l 


5.0542 


4435 


= h-h 


25.9961 


280 




60.4039 


165 


=h+h+h 


5.3804 


119 


= /i4 - h 


26.2959 


441 


= /l6 


63.7054 


152 


= 2/2 + h 


5.5458 


118 


= 2/2 -h-h 


26.4532 


524 


= 2/2 - /l 


63.7415 


126 


= 2/2 + h 


6.2363 


139 


= 3/i - 2/2 


26.5422 


148 




64.1969 


292 


= 3/2 


6.3155 


678 


= /5 - /l 


27.2590 


105 


= /4 + /5 - /l 


65.3790 


628 


= 4/1 


6.7278 


108 


= 3/l - /2 - h 


27.6352 


481 


= 3/i - /2 


69.9416 


418 


= 3/1 + h 


9.1253 


116 


= 2/3 - 2/1 


27.7144 


172 


= /2 + /5 - /l 


69.9778 


241 


= 3/1 + /4 


9.1615 


112 


= /3 + /4 - 2/i 


28.0907 


238 


= 3/l - /4 


70.4332 


318 


= 3/1 + h 


9.6169 


204 


= /2 + h ' 2/1 


28.1268 


426 


= 3/1 - /3 


71.6944 


70 


= 3fi+h 


10.0293 


361 


= 2/i - h 


29.6549 


393 




74.5043 


94 


= 2/1 + 2/3 


10.1085 


233 


= 2/2 - 2/1 


29.6621 


158 




74.5404 


132 


= 2/1 + /3 + /4 


11.2905 


2066 


= 2/i - /2 


30.5486 


112 


= 2/6 


74.9959 


201 


= 2/l + /2 + /3 


11.7459 


1068 


= 2/1 - /4 


31.6190 


196 


= /l+/6 


75.0320 


79 


= 2/l + /2 + /4 


11.7821 


1844 


= 2/i - h 


32.1979 


141 


= 2/l + /3 - /2 


75.4874 


271 


= 2/l + 2/2 


12.0343 


41 


= 2/1 - /12 


32.6895 


12311 


= 2/1 


76.7487 


67 


= 2/1 + /2 + /5 


12.2342 


62 


= 2/i - /ii 


32.7252 


112 


= 2/i +h~ h 


80.0501 


99 


= 2/2 + /l + /3 


12.4471 


121 


= 2/1 - /lo 


33.1811 


148 


= 2/1 + /2 - /3 


80.0863 


110 


= /l + 2/2 + /4 


12.6709 


107 


= 2/i - h 


33.4618 


272 


= h+h 


80.5417 


142 


= 3/2 + /l 


12.9060 


81 


= 2/1 - /8 


36.1282 


175 


= h+h 


81.7237 


182 


= 5/1 


14.5919 


108 


= /l+/3-/5 


36.3633 


216 


= h+h 


81.8029 


41 


= /l + 2/2 + /5 


15.0835 


262 


= /l+/2-/5 


36.5871 


270 


= fl+ /lO 


85.5959 


68 


= 4/2 


15.2743 


692 


= /6 


36.7935 


112 


= /l + /l3 


86.2864 


146 


= 4/1 + /3 


15.8532 


425 


= /l + /3 - h 


36.8000 


136 


= /l + /ll 


86.3225 


77 


= 4/1 + /4 


15.8893 


180 


= /l+/4-/2 


37.0000 


74 


= /l + /l2 


86.7780 


134 


= 4/1 + /2 


16.3086 


230 


= /l+/3-/4 


37.2410 


192 




88.0392 


29 


= 4/l+/5 


16.3618 


155 




37.2521 


4143 


= fi+h 


90.8490 


46 


= 3/i + 2/3 


16.3807 


285 


= /l + /4 - /3 


37.2883 


3158 


= fi+h 


90.8852 


83 


= 3/l + /3 + /4 


16.7999 


141 


= h+h-fA 


37.7437 


5730 


= /l+/2 


91.3406 


120 


= 3/l + /2 + /3 


16.8363 


369 


=h+h-h 


37.9345 


183 


= /5+/6 


91.3768 


32 


= 3/1 + /2 + /4 


17.1170 


871 


= h 


38.0698 


170 


= /l + /l4 


91.8322 


97 


= 3/i + 2/2 


17.6060 


212 


=h+h-h 


39.0049 


795 


= /l + /5 


93.0934 


38 


= 3/l + /2 + /s 


19.3224 


424 




40.1192 


146 


= /l + /l5 


95.9033 


32 


= 2/l + /2 + 2/3 


19.7078 


341 




40.6909 


295 


= /3+/8 


95.9394 


28 


= 2/l + /2 + /3 + /4 


19.7835 


636 


= h* 


41.4176 


128 


= /2 + /9 


96.3949 


91 


= 2/2 + 2/i + /3 


19.8051 


152 




41.8148 


494 


= 2/3 


96.4310 


42 


= 2/1 + 2/2 + U 


20.0186 


814 


= h* 


41.8509 


760 


= /4 + /3 


96.8864 


68 


= 3/2 + 2/i 


20.2424 


959 


= /lO* 


42.3064 


1139 


= /2 + /3 


98.0685 


57 


= 6/1 


20.2487 


126 




42.3427 


233 


= /2 + /4 


98.1476 


27 


= 2/1 + 2/2 + /5 


20.4487 


508 


= /l3 


42.6408 


103 


= /l + /l6 
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Table 5. Properties of the best TDC model, obtained by perturbing the spectroscopic parameters and convection efficiency. 



Mass 
Mq 


Radius 

R© 


K 


Age 

Myr 


log 9 

(cgs) 


log(L/LQ) 


X 


Z 


"MLT 


UiQV 


1.56 


2.425 


7110 


1829 


3.862 


1.130 


0.7465 


0.0071 


2.0 


0.2 



Table 6. Characteristics of the best representative rotating model. The first six columns have their usual meanings, then headings are: 
mean density, equatorial rotation velocity, the ratio of the equatorial rotation velocity to the keplerian orbital velocity at the stellar 
surface, inclination of the stellar rotation axis to the line of sight, fundamental radial mode frequency f(F0), first overtone radial mode 
frequency f(Fl), the base-10 logarithm of the period of FO, and the period ratio of FO and Fl. The model has metallicity [Fc/H]=— 0.52, 
convective efficiency, ajviLT = 0.5, and overshooting dov = 0.3. Sec ^ l5.2l for model selection criteria. 



Mass 


Radius 

R© 


Tcff 
K 


Age 
Myr 


logs 

(cgs) 


P 


kms^-"- 


i 
dcg 


f(FO) 
d-i 


f(Fl) 
d-i 


log P{FO) 


P(F0)/P(F1) 


1.53 


1.764 


7256 


1465.93 


4.129 


0.392 


34.18 


47til 


16.327 


21.1 


-1.213 


0.773 
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